Correlation effects during liquid infiltration into hydrophobic nanoporous mediums 
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Abstract 



Correlation effects arising during liquid infiltration into hydrophobic porous medium are 
considered. On the basis of these effects a mechanism of energy absorption at filling porous 
medium by nonwetting liquid is suggested. In accordance with this mechanism, the absorption 
of mechanical energy is a result expenditure of energy for the formation of menisci in the 
pores on the shell of the infinite cluster and expenditure of energy for the formation of liquid- 
porous medium interface in the pores belonging to the infinite cluster of filled pores. It 
was found that in dependences on the porosity and, consequently, in dependences on the 
number of filled pores neighbors, the thermal effect of filling can be either positive or negative 
and the cycle of infiltration-defiltration can be closed with full outflow of liquid. It can 
occur under certain relation between percolation properties of porous medium and the energy 
characteristics of the liquid-porous medium interface and the liquid-gas interface. It is shown 
that a consecutive account of these correlation effects and percolation properties of the pores 
space during infiltration allow to describe all experimental data under discussion. 

1 INTRODUCTION 

Energetics of "nanoporous medium-nonwetting liquid" systems is one of the new directions 
in basic and applied research (see e.g. [THE])- In the simple model of a porous medium in the 
■<sj- ! form of cylindrical channels, threshold pressure is described by the Laplace- Washburn equation 
p = 2a/ R\cosa\, where a is the surface energy of the liquid, R is the pore radius, and a is the 
contact angle (for a nonwetting liquid, a > 90°). For filling nanometer-sized pores by a nonwetting 
liquid with a surface energy of 0.05 0.5 J/m 2 , this threshold pressure is 10 2 -j- 10 3 atm. When the 
liquid passes from the bulk state to a dispersed state in pores of the nanoporous medium with a 
specific volume of ~ 1 cm 3 /g, the absorbed and accumulated (returned when the liquid flows out) 
| energy can reach 10 -r- 100 kJ/kg. This value is one order of magnitude higher than the energy 
observed for widely used materials, such as polymer composites and alloys with the shape memory 
effect [10HT2]. 

Among the systems under investigation are silochromes, zeolites with liquid metals, hydropho- 
bized silica gels, and zeolites with water and aqueous solutions of organic compounds and salts. In 
recent years, hydrophobized nanoporous media have become available owing to the development 
of the method used for modifying the surface of nanoporous media, generally, with alkyl chlorosi- 
lanes [3ll6ll71 fT2Tl20] . To date, nanoporous media with different pore shapes, porosities, specific 
surface areas, specific volumes, average pore radii, and pore size distributions have been stud- 
ied [THTl [9rH~2 | ITTrl50] . The investigations performed thus far have been concerned primarily with 
equilibrium properties. Experiments have been carried out at a low compression rate of the system 
when the rate of increase in the pressure (p) in the liquid-porous medium system is (10 -3 — 1) 
atm/s. In the infiltration-defiltration cycle, there is a hysteresis, so that the threshold pressure of 
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infiltration is higher than the pressure of defiltration. Moreover, the majority of the systemsare 
characterized by the phenomenon of nonoutflow of a nonwetting liquid when a part of this liquid 
remains in the porous medium as the excess pressure decreases to zero. The absorbed energy is 
determined by the product of the volume of filled pores and the difference between the infiltration 
and defiltration pressures. In frameworks the model of cylindrical channels, these pressures are 
described by the Laplace- Washburn equation with different contact angles. 

The revealed difference between the infiltration and defiltration pressures and the absorption of 
the mechanical energy observed in the infiltration-defiltration cycle due to the pressure hysteresis, 
as a rule, have been explained by the hysteresis of the contact angle; however, the mechanism 
responsible for the appearance of the latter hysteresis has remained unclear [4ll6ll7l lT9l - [2~2] . 

It has been established that the infiltration and defiltration pressures depend on the tempera- 
ture and that, for the porous medium with a disordered structure of pores, the defiltration pressure 
increases (by several factors) with an increase in the temperature from 280 to 350°K, whereas the 
infiltration pressure decreases only slightly (by less than 10%) or remains constant [2~T | |2~2" } |2~T] . This 
means that, during infiltration and defiltration, the phenomenological contact angles differently 
depend on the temperature. For zeolites, the revealed temperature dependences exhibit a more 
complex behavior; moreover, the volume (V) memory effect can be observed with an increase in 
the temperature and its subsequent decrease. It is worth noting that this effect is one order of 
magnitude (in AV/V) stronger than that observed for known alloys and composites [TO ] ITT ] . 

At present, there exist several hypothesis regarding the nature of the contact angle hysteresis. 
This hysteresis has been attributed to the rough surface of pores, the chemical inhomogeneity of 
the surface, and the dependence on the direction of the liquid motion [T91I20] . 

However, these conceptions cannot explain the different dependences on the temperature of 
infiltration and defiltration pressures, and hence the temperature dependence of contact angle. 

In the framework of the concept that the porous medium is a system of cylindrical channels, 
the absorbed energy is expended for forming a liquid-porous medium surface, which appears in 
the course of infiltration and disappears during defiltration at different pressures due to different 
contact angles (21 El SSI ST] - In case of the closed hysteresis loop in the infiltration - full- 

defiltration process initial and finite states of the system are similar, the internal energy change 
AEdcie = § dE = and the work done to perform the filling of the porous medium should be equal 
to the thermal effect § dA = § dQ. The measurements carried out in [T91I5T] showed that rise of 
temperature in the hydrophobic silica gel-water systems under investigation during infiltration- 
defiltration cycle was <10~ 3 K. On the other hand, it was found in |13] that when one-third of the 
porous medium volume of a similar system was filled and the liquid-porous medium interface area 
change was not equal to zero, the temperature did not increase within the limits of error (< 0.1K), 
while the estimation provides the temperature increase by AT = 0.8K. 

In [2J the dependence of the thermal effect AQ, which is accompanied by the thermal effect 
due to elastic compression of water and the porous medium, on the filling degree in the hydropho- 
bic silica gel-water system was measured. The author found that during infiltration this porous 
medium with ip m 0.22, the heat generation took place (AQ < 0), whereas ^ < for water [51] 
and, consequently, in case of independence of the contact angle from temperature during the for- 
mation of liquid-solid interface with area AS one would expect the AQ = —j^AS value to be 
positive, i.e. heat absorption shall be observed. In [2J, it was experimentally found that, as the 
degree of infiltration increases, the quantity \AQ\ reaches a maximum at 9 « 0.6 and decreases 
to zero after the complete infiltration of the porous medium: AQ = 0. This result led the author 
of [52] actually to the conclusion that the first law of thermodynamics is violated in the process 
under investigation. 
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Thermal effect associated with the infiltration of the modified porous medium was also observed 
|47ll4"8] . It was found that for the investigated porous media with (p = 0.4 and (p = 0.6 the value 
AQ < reduces with increasing degree of filling, reaching a minimum at complete filling. 

Note that the energy absorption can not be to all appearance explained by viscous dissipation, 
since, as it ascertained in [59], change in viscosity of nonwetting liquid (aqueous solution of CaCLj) 
by 7 times does not alter the dependence of the filled volume on time, the threshold infiltration 
pressure and rate of filling hydrophobic silicagel. 

Statistical theory of fluctuations allows us to generalize the Laplace- Washburn equation in case 
of filling the pores in the system of interconnected pores [HE]. The proposed approach makes it 
possible to analyse the probability of the fluctuation infiltration pores with the size R near the 
threshold pressure. If the probability of w ~ exp(— SA/T) increases continuously with growing 
pressure, then the pore can be filled with liquid (5 A is work expended for infiltration of one pore, 
5A(p, a, 5a, R)). The value 5Ai n takes into account a contribution of formation a liquid-porous 
medium interface and a surface of menisci at the mouths of the filled pores. Threshold infiltration 
pressure is determined from the condition 5Ai n (p out , a, 5a, R) = 0, and the defiltration pressure is 
determined from the condition 5A out (p out , a, 5a, R) = 0. This makes it possible to calculate contact 
angles during infiltration and defiltration, to describe the dependences of the volume of system 
liquid-porous medium on the pressure and calculate the volume of confined liquid after defiltration. 

In [53j the energy absorption during the infiltration-defiltration process within the framework of 
fluctuation model is associated with the energy of formation of menisci. However, in this model is 
assumed, that the work of filling the pores (liquid defiltration) and, consequently, the probability of 
infiltration (defiltration) does not depend on the degree of filling of the porous medium. Therefore, 
if during infiltration and defiltration the contact angle is independent of the filling degree, the area 
of menisci which are formed should be equal to the area of menisci disappeared, and therefore the 
formation energy at full infiltration (defiltration) should be equal to zero. It also does not allow 
to explain the results of [2||4^||4T]. 

Thus, the experimental data currently known seems to be self-contradictory and well-known 
traditional mechanisms of energy absorption do not make it possible to explain the infiltration- 
defiltration hysteresis, temperature dependences of the infiltration and defiltration pressures and 
special features of heat generation during the filling of nanoporous media with nonwetting liquid 
which are observed. 

As shown in this paper (Section 2), to describe the discussed phenomena we should take into 
account the structure of pores in disordered porous medium. 

During infiltration into porous medium with an increase in the number of filled pores, sur- 
rounded by empty filled pores a formation meniscus in pore mouths is changed by theirs disap- 
pearance on the average ensemble of pores with increasing number of full filled pores surrounded 
by empty loadable pores. In this case expenditure of work on formation meniscus are replaced in 
the gain energy during theirs disappearance. Thus, the correlation effect of the relative position 
of filled and empty pores must be taken into account at describing the fluctuation infiltration. 
In addition, the porosity of the medium determines the possible neighbors number of pores and 
possible macroscopic set of filled pores in the percolation cluster, through which liquid flow from 
the surface of the porous medium to fill the pores is possible. This means necessity of taking into 
account another correlation effect - the spatial arrangement of pores in the medium, arising as a 
result of probabilistic system realization of interconnected pores in the infinite percolation cluster. 

Consequently, it is possible to sugges the mechanism of energy absorbtion during filling of 
the nanoporous medium with nonwetting liquid is proposed (Sections 3, 4). In accordance with 
this mechanism, the mechanical energy absorption is a result of expenditure of energy for the 
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formation of menisci in the pores on the shell of the infinite cluster and expenditure of energy 
for the formation of liquid-porous medium interface in the pores belonging to the infinite cluster 
of filled pores. It has been demonstrated (Section 4), that the infiltration-defiltration cycle can 
be closed with the complete defiltration of the liquid and the reproducibility of the cycle when 
the specific relationship between the percolation properties of the porous medium and the energy 
characteristics of the liquid-porous medium and liquid-gas interfaces is satisfied. It turned out 
that, depending on porosity and, consequently, the number of nearest neighbors of the filled pores, 
the thermal effect can be either positive or negative. This makes it possible to explain the known 
experimental data [4"3~] , |4"6 | |4"T ] mentioned above and to describe the temperature dependences of 
the infiltration and defiltration pressures [21~ ll22ll24] . The description of experimental data on the 
basis of the proposed approach is given in sections 5, 6. 

2 THE MODEL OF A POROUS MEDIUM. INFILTRA- 
TION FLUCTUATIONS 

Let us consider a disordered porous medium infiltrated with a nonwetting liquid. It is assumed 
that the half-width of the pore size distribution 5R satisfies the inequality SR/R < 3, so that the 
fulfillment of this inequality ensures the independence of the percolation threshold from the radii 
of pores [54j. 

Filling of the porous medium is a process of liquid infiltration into the disordered porous 
medium which contains pores of different size. It is assumed that the size of the porous medium 
a is much more than the maximum size of the pores -R max ( a > 10 3 -R m ax |55j) so that the porous 
medium can be regarded as infinite. Obviously, the infiltration of all the pores can take place only 
when the pores are connected with the surface and form a connected system. Consequently, the 
filling of the porous medium can take place only when the pore system in it is far beyond the 
percolation threshold (<p > ip c ), ip is its porosity equal to the ratio of pores to the porous medium 
volume, tp c is the percolation threshold, which is the characteristic of the porous medium. For 3D 
systems, the percolation threshold p c = 0.18 [56j. At the same time, the connectivity of pores with 
one another is the result of the formation of infinitely large clusters of pores at ip = <p c . Figure [1] 
shows the dependence of the probability normalized to unity of a pore belonging to the infinitely 
large cluster on porosity ip [56,57j. 

It can be seen from Fig. [I] that, in the vicinity of the percolation threshold ip c , only a small 
number of pores (~ 1%) belong to the infinite cluster; therefore, in this case, only a small fraction 
of these pores, as well as pores that belong to the finite clusters connected with the boundary of 
the porous medium, can be infiltrated. At increasing porosity and for (p ^> tp c P(<p) — > 1 and, 
consequently, the pore space becomes homogeneous due to the growth of the infinitely large cluster 
of pores. Under these conditions, the infiltration of the porous medium can be described as the 
infiltration of an infinite cluster of pores. It is this infiltration that will be considered in the present 
paper. 

It is assumed that thermal fluctuations at the pressure p in the vicinity of the infiltration 
threshold of the porous medium bring about the formation of macroscopically small regions in the 
form of clusters consisting of N pores filled with a liquid. Each cluster arises at the boundary of 
the porous medium and, in the view of the boundedness of the pore volume, grows through the 
attachment of the other filled pores to it. We believe that, at the beginning of the growth, each 
cluster can be considered a system of branched chains consisting of filled pores. In the course 
of infiltration of the porous medium with the liquid, the external pressure does the work. This 
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Figure 1: The dependence of the probability normalized to unity of a pore belonging to the 
infinitely large cluster on porosity ip 

process is accompanied by the formation of energetically unfavorable surfaces of both the menisci 
of the liquid in pores and the liquid-porous medium interfaces. Moreover, the state of the gas 
in the pores and the elastic state of the porous medium change as well. If the adiabatic work of 
formation of an infiltration fluctuation is 5A(N) and the energy of dissipation due to the friction 
can be disregarded (see [55J), the probability of the formation of a fluctuation can be written as 
w ~ woexp(5S) [58], where 5S = —5A/T is the fluctuation of the entropy. Therefore, an increase 
in the quantity 8 A (5 A > 0) with an increase in the number of pores N in the cluster leads to 
a decrease in the fluctuation probability. This corresponds to the thermodynamic stability of the 
initial state of the system. For 5A ~ T, the infiltration fluctuation can increase. In this case, the 
system becomes unstable and the liquid begins to infiltrate the porous medium. 

The infiltration of a nonwetting liquid under the pressure pi n a porous medium requires a 
certain amount of work to be done to fill the pores of the porous medium. For this purpose, it 
is necessary to overcome a certain critical pressure, which is the Laplace pressure p c (R) ~ ^ for 
an isolated pore with characteristic size R. 5a = a sl — a sg , where a sl is the surface energy of the 
solid-liquid interface, <r lg = a is the surface energy of the solid- gas interface, 5a = a sl |cosa|, a is 
the contact angle. 

At a pressure lower than the critical value pi n , the adiabatic work satisfies the inequality 
5A(N) > at any value of iV and the fluctuation probability decreases with an increase in N. 
Therefore, the fluctuation probability is equal to zero for any macroscopically large number of 
pores. Fluctuations of finite length arise, but no infiltration of the porous medium occurs. At a 
pressure in the vicinity of the critical value p in , the work is 5 A ~ T and thermal fluctuations in the 
system can lead to the formation of clusters from N pores. At a pressure p > pi n , the infiltration 
of individual pores becomes energetically favorable because the quantity 5A(N) is negative. Since 
the work is 5 A ~ N, the fluctuation probability at a pressure p > pi n is i« ~ 1. The pressure 
difference p — Pi n causes the liquid to move in the porous medium. 

Now, we consider a porous medium immersed in a nonwetting liquid under an external pressure 



p, which does the work in the course of infiltration of the porous medium. Let 5 A (p) be the work 
expended for providing the fluctuation infiltration of one pore. According to [5], the expression for 
the work 5 A (p) with due regard for the formation of menisci can be written in the form: 

5A(p) = -pV + aS rn + 5a(S-S m ). (1) 

Here, V is the volume of the pore, S is the surface area of the pore, S m is the surface area of the 
menisci, and a is the surface energy of the liquid. 

For a spherical pore with the radius R, the work 5 A (p) can be represented in the form [5]: 

A , r,\ 47rR3 

5A = A(p, R — — , 

3 (2) 

35a - 
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where rj is the ratio of the meniscus surface area to the pore surface area. 

A similar expression can be written for a cylindrical pore. For example, we write the following 
expression for the work expended for providing the fluctuation emptying of a cylindrical pore by 
nucleation with the nucleus length L and the radius R: 



5A = A (p, R, L) A/3ttR 



•3 



A(p,R,L)=p{ 



3x _ , 35a 



a 



(3) 



2 ' R 
where x = L/2R. 

It follows from relationship fl2]) that the infiltration of the pore becomes energetically favorable 
at a pressure higher than the critical value pi n . The critical pressure is determined by the equality 
of the work on the fluctuation infiltration of the pore to zero. An analysis of the critical pressure 
for a spherical pore was performed in |5]. Expression fl3]) implies that, devastation of the pores 
becomes energetically favorable at a pressure less than the critical pressure p ou t, which depends 
on the nucleus length L and the nucleus radius R and, at L — > oo, transforms into the Laplace 
pressure p ou t ~ This means that, in the work expended for infiltrating the pore, the change 
in the surface energy of the pore dominates over the energy expended for forming menisci. 

For characteristic values of the parameters of the porous medium and the liquid a ~ 5a ~ 0.5 
kJ/m 2 , R ~ 10 nm, and, \ ~ 10~ 4 atm -1 , when the pressure deviates from the value pj by 5p = 
— 10~ 2 pi, the work changes from 5 A ~ T to ~ 1 eV. In this case, the change in the probability of 
infiltration fluctuation at the temperature T = 400K is equal to ten orders of magnitude. Therefore, 
for systems with the aforementioned characteristic parameters, the probability w changes abruptly 
from to 1 in a narrow pressure range (5p/p ~ 10~ 2 ) in the vicinity of the pressure pi n . It should 
be noted that the inclusion of the gas filling the pores in the analysis leads to the appearance 
of an additional contribution to the work 5A g . The value of this contribution 5A g during filling 
one pore with the volume V can be estimated under the assumption that it is equal in the order 
of magnitude to the specific evaporation energy of the liquid fio multiplied by the density of the 
gas p g at the pressure p, i.e., 5A g ~ fi p g V. For water, we can write p ~ 2.2 • 10 3 J/kg and 
Po ~ 1-0 kg/m 3 . Assume, that the filling of the porous medium occurs at a pressure of the liquid 
Pi n ~ 1 • 10 7 J/m 3 . In this case, the additional contribution to the work due to the presence of 
the gas in the pores is given by 5A g ~ PoP g V ~ 10~ 2 • 5 A. Nominally, accounting of the influence 
of gas reduces to the appearance of additional summands in equations ([2]) and ([3]), each of them 
proportional to the gas pressure at a given temperature. The value of these summands A g i n and 
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Agout for typical infiltration pressure pi n ~ 200 atm and p on j ~ 1 — 10 atm, for pore radius R « 10 
nm is A S j n ~ A g ~ v4 sout ~ 10 _2 eV in order of magnitude. Therefore, the influence of the gas 
in the pores on the infiltration of the porous medium can be ignored at considering the filling 
of the pore, but accounting of these influence may be important at considering the effects of the 
defiltration liquid the cylindrical pores. 

Since the pore can be either filled (probability w ~ 1, 8 A (p) < 0) or empty (probability w = 0, 
S A (p) > 0), the normalized probability can be written in the form [5j 

w l (p) = [l+exp(5A/T)]- 1 . (4) 

Note, that relations ([3]) and (j3J) explain the obtained in [6|,[7] experimental data for the depen- 
dence of infiltration pressure and defiltration pressure on the pore size. 

Pores are not isolated in a porous medium, but they are connected with one another by throats 
(mouths), in which menisci are formed during infiltration of a certain pore. Thus we can assume 
that the medium subjected to infiltration is the heterogeneous medium which consists of full and 
empty pores playing the role of white and black spheres, respectively, in the percolation theory |56| . 
Such medium can experience percolation transition occurring via the formation of clusters of 
accessible pores and followed by infiltration of nonwetting liquid into such formations. [25] • In 
addition to that, filling of the macroscopic volume of a porous medium occurs by infiltration in 
the infinitely large cluster of accessible pores [55J . 

Below, we will consider the infiltration of pores located on the shell of an infinite cluster 
consisting of filled pores. In this case, the condition 5A = determines the pressure necessary for 
the infiltration of a pore on the shell of the infinite cluster of filled pores. 

3 WORK AND THERMAL EFFECT IN THE INFILTRATION- 
DEFILTRATION CYCLE 

Let us calculate the work 5 A and the thermal effect AQ for an arbitrary degree of infiltration 9 
of the porous medium. The thermal effect AQ in filling of a porous medium by nonwetting liquid 
comprises the thermal effect AQ p due to the formation of the liquid-solid interface , the thermal 
effect AQ W related to formation-disappearance of menisci and the thermal effect AQ U related to 
the compressibility of the nonwetting liquid-nanoporous medium system. 

AQ = AQ P + AQ W + AQ U . (5) 

The AQ P , AQ W and AQ U values can be calculated using thermodynamic relations [58J, which 
determine the thermal effect in formation of the surface AQ S : 

AQ S = -^AS. (6) 
Here, AS is the change in the system surface. 

To calculate the thermal effect AQ S let us suppose that each pore in a porous medium has z 
nearest neighbours and pores contact each other by throats, each of which has an area S z . If an 
empty pore contacts a full one, the meniscus is formed in the throat. A pore in a porous medium 
can be filled only if liquid can reach it. In compliance with above-mentioned assumptions this 
condition can be satisfied by formation of an infinitely large cluster of filled pores. In this case, 
only those pores which belong to the shell of the infinite cluster will be filled. It is possible to 
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show that the contribution of filled finite-size clusters (which liquid can reach via filled clusters 
contacting with the surface of the porous medium)to the filled volume is small. Distribution f(N) 
of the number of pores in clusters of finite size near the percolation threshold is determined by 
the scaling dependence f(N) ~ , r = 2.2 [56]. Therefore the bulk of the cluster contains one or 
more pores, which are mostly not associated with the surface of the porous medium. Liquid cannot 
reach such pores and, consequently, they are not filled at 9 ~ 9 C . Taking it into consideration, we 
can represent the thermal effect 5Q related to filling of one pore as: 

6Q = ~ T ^ {S ~ zSz) ~ T ^ zS * w W> ( ? ) 

S = 4nR 2 is the area of the surface of a pore with radius R, z is a number of nearest neighbours, 
W(9) is the difference (averaged over the ensemble of pores) between the numbers of menisci before 
and after the infiltration of the pore per the nearest neighbor of the infinite cluster. 

Considering that filling of a porous medium is the result of formation of the infinitely large 
cluster of filled pores and taking into account the normalized probability P(9), we obtain that the 
quantity of heat per a pore released in the process of filling the porous medium to the degree of 
filling 9 can be written as: 

6 

AQ P (9) = -T-£ J < {S - zS z ) > 9P(9)d9, 
o 

oo 

<S-zS z >= J dR{S - zS z )f(R), (8) 

Q 

6 

AQ W (9) = -T^ J < zS z W(9) > d9. 
o 

Here, f(R) is the function normalized to unity of the size distribution for pores. For a disordered 
porous medium average values in (JHJ) can be calculated in the framework of a specific model of a 
porous medium. We will use the model of randomly arranged spheres in which pores represent 
randomly arranged spherical holes [59] ■ This model does not take into account correlations in 
location of pores with different radii in accordance with assumption made about the narrowness of 
size distribution for pores (5R) <C R- In accordance with the model, the average number of nearest 
neighbours z, associated with porosity of medium <p, and the area of a throat can be written in 
forms [59] : 

z = z = — 8 ln(l — ip), 

s. = (9) 

256 

Using expression (JHJ) and taking into account that (SR) R, we obtain from flU]): 



aq p (9) = -t^(i - v)^R 2 J ep(e)de, 



e 

da -o f (10) 
AQ W = -T V — 4vri? 2 y 

o 

9tt 2 

ri = z . 

' 1024 

To calculate W{9) we consider an empty pore located on the perimeter of the infinite cluster 
of filled pores. Let us suppose that this pore contacts the infinitely large cluster of filled pores via 
n throats. Thus, menisci are formed in all n mouths and menisci are absent in the remaining z — n 
throats. After filling this pore, menisci which were there at the beginning of infiltration disappear 
and the number of menisci will be equal to z — n. In this case we can write W{9) as: 

W{9) = £ (P(0)r(i - ey-^ z ~ 2n zl (ii) 

z — ' z n\(z — n)\ 

71=1 v ' 

The first factor under the summation sign determines the probability that an empty pore 
contacts the infinite cluster of filled pores n times, the second factor is the probability of finding 
the empty pore close to the infinite cluster, provided that this pore is surrounded by z — n empty 
pores and therefore has z — n throats. The third factor determines the difference between a relative 
number of menisci after (z — n) and before (n) the filling of the pore. The binomial coefficient takes 
into account variants of allocation of n menisci on number of pore nearest neighbours. Note that 
the obtained expression coincides with the full perimeter of the infinitely large cluster calculated 
in |60j, if the third factor is substituted for unity. 

The sum in Eq. (TTTj) can be calculated analytically: 

W{9) = (6 2 - 26 - P(6) + 1 + 6P{0)){P{6) -6 + I) 2 " 1 - (1 - 6) z+l . (12) 

Relations (Q, (Tl0|) . ( ITS]) determine the thermal effect during filling a porous medium with 
porosity p to the degree of filling 6. 

The work AA expended for infiltrating the porous medium to the volume determined by the 
fraction 9 and the corresponding work expended for infiltrating one pore 8Ai n (9) can be calculated 
from the thermodynamic relationship A = J adS [58]. By using expression (j2J) for spherical pores, 
we obtain: 

5A in (6,p) = -P^z- + 5a(l - t])AtiR 2 + ar)AnR 2 W{e), 

A A = AA P + AA W , 
e 

AA P {6) = 5a (1 - ri)47iR 2 J 6P(6)d6, (13) 

o 

e 

AA W = ariAnR 2 J W{6)d6 
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The sum of heat ( 1101) and work ( 1131) determines the change in the energy of isothermal infil- 
tration of the porous medium: 



AE = AE P + AE W , 



AE P = (5a - r^)(l - vHnR 2 [ 9P(9)d9, 

dl i (14) 

e 

AE W = (a- T^) V 4nR 2 J W(9)d9. 

o 

It follows from (THj) that the change in energy of the system during filling the porous medium is 
determinded by the specific surface energies a and 5a, geometric properties of the porous medium 
and the evolution of the infinite cluster of filled pores, which depends on the properties the disor- 
dered porous medium (Fig. [T]). 

For the calculation of the work 5A{9) and the thermal effect AQ V arising upon the defiltration 
of the liquid from the porous medium, it should be noted that, in the infiltrated porous medium, 
the defiltration of the liquid leads to the formation of empty pores surrounded by at least one 
filled pore connected through other filled pores with the surface of the porous medium. As in 
case of filling, the formation of an empty pore goes with the change in surface energy of the 
liquid-solid and liquid-gas interfaces as well. The change in surface energy is associated with 
formation-disappearance of menisci [5j. 

Taking into account this fact, the work expended for emptying the pore in the porous medium 
5A(8) with the degree of infiltration 9, the work for defiltration liquid AA V from 1 to 9 the degree 
of infiltration, and the thermal effect 5Q V associated with the defiltration of the liquid from one 
pore in the porous medium can be written in the form: 

4-7T 3 

5A 0Ut {9) = p Sa{l - ri)4nR 2 + ar}AirR 2 W x {9), 



(15) 





AA v p (9) = Sa(l - r])AnR 2 J 9' d& , 

o 

9 

AA V W = -ariATrR 2 J W 1 (9 , )dB', 
o 

SQ V = ~T^(S - zS z ) + T^zS z Wi(9). 

Relations (J7]) and (lT5"j) differ in sign of the last term and the functions W(9) and W\{9) which 
determine the difference per one nearest neighbour between the number of menisci before and 
after infiltration (defiltration) in pores. 

In contrast to the case of infiltration, the defiltration of the liquid occurs initially through the 
formation of individual empty pores and clusters of empty pores with a decrease in the pressure 
and, after the infinite cluster of empty pores is formed, through the formation of pores on the 
shell of this cluster. Upon the defiltration when the low degree of infiltration is reached, the liquid 
can be retained in the porous medium if it is contained in individual pores or clusters of filled 
pores surrounded by empty pores with smaller sizes from which the liquid defiltrated at higher 
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pressures. However, as the number of neighbors of empty pores increases, according to relationships 
( I15jl (see below Figs. [JJ|3]), the defiltration of the liquid becomes energetically more favorable; i.e., 
it should proceed at higher pressures. Therefore, the quantity W\(9) should be defined as the 
difference (averaged over the ensemble of pores) between the numbers of menisci before and after 
the emptying of the pore on the shell of the system of empty pores. Taking it into consideration, 
calculation of Wi(9) gives: 

Wi{6) = (29 2 - 30 + 1) - (1 - 6) z+1 . (16) 

The thermal effect, the work, and the change in the energy during the defiltration, when the degree 
of infiltration varies from 1 to 9, can be written in the form similar to relationships (fTUj) . (fT3|) . and 



aq;(9) = -T^ii - v)^R 2 J ode, 

o 

e 

AQl = Tv^rcR 2 J W 1 (9)d9, 
o 

6 

AA v p {9) = 5a{l - t])AttR 2 J 9d9, 

o 

6 

AA V W = -ariATcR 2 J Wx{6)d9, 
o 

e 

AE v p = (6a - T^)(l - V )47rR 2 J 9d9, 

o 

e 

AEl = -(a - T^v^R 2 J W x {e)d9. 



(17) 



o 

Expressions (1T7|) . like relationships (TTDj) . (lT3"j) . and (I14p . are valid for the case of an isothermal 
process. This implies that they can be used for describing experiments in the case where the 
characteristic time of heat transfer (removal) tq is considerably shorter than the characteristic 
time Ty of change in the volume of the nanoporous medium-nonwetting liquid system: tq <C Ty. 
At tq > Ty, the temperature and, correspondingly, the quantities a, 6a, da/dT, and d6a/dT 
become dependent on the time and, hence, on the degree of infiltration 9. In this case, they should 
be introduced under the integral sign in relationships (TL41) and (|T7|) . The inequality tq Ty impose 
constraints on on velocity of compression system under the study of its equilibrium properties. 



4 CONDITIONS FOR THE CLOSED CYCLE 

It can be seen from relationships ( JT4|) and ( ITT]) that, if after the increase in pressure and 
the infiltration of all pores of the nanoporous medium with a liquid and the subsequent decrease 
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in pressure and the complete defiltration, the system reverts to its original state, the following 
relationship should hold true: 

1 1 
(5a - T^)(l -n)J (6P(6) - 9)d9 + (a - T^) v J {W 1 {9) + W{9))d9 = AE, 

o o 
AE = 0. 

Here, AE is the change in the internal energy of the system upon the transition from the initial 
state to the final state in the course of infiltration and defiltration. 

Expression (1TB]) relates the energy parameters of the liquid-solid and liquid-gas interfaces and 
the macroscopic characteristics of the porous medium, such as the porosity ip, the structure of 
the percolation cluster P(9), and the quantities W and W\, which determine the dependence on 
9 for the surface of menisci at the mouths of filled pores on the shell of an infinite cluster of 
filled pores in the case of infiltration and for the surface of menisci in pores on the shell of all 
clusters of empty pores (including the infinite cluster) in the case of defiltration. It follows from 
relationships (118j) that the absorbed energy in contraction and expansion of the system in the 
closed cycle is dependent on the quantities a and 5a and equal to the total heat released upon the 
formation and disappearance of the liquid-solid and liquid-gas surfaces. This heat is determined 
by the independent quantities, namely, the derivatives of the quantities a and 5a with respect to 
temperature. The integrals in expression (fT8l) account for the different paths of the system in the 
course of infiltration and defiltration. In the closed cycle, according to relationships (11 Op . (I13p . 
( TTTjl . and ( TT8l) . during the infiltration and defiltration, the system undergoes different sequences of 
equilibrium states that differ in macroscopic sets of filled pores. 

In particular, the infiltration of the porous medium according to expression (TT8]) is accompanied 
by an increase in the number of pores that belong to the infinite cluster of filled pores and by a 
change in the number of menisci in pores on the shell of this cluster. The defiltration of the porous 
medium is accompanied by an increase in the number of pores in all clusters (including the infinite 
cluster of empty pores and single pores) and by a change in the number of menisci on the shell of 
the entire system of empty pores. As follows from relationships (Q, (fTUj) . and (fT6|) . these sequences 
of states depend on the porosity <p and the number of neighbors z in the system of connected pores. 
Therefore, in terms of the percolation theory and the model under consideration, the contact angle 
hysteresis is associated with different (in infiltration and defiltration) macroscopic properties of 
systems of filled and emptied pores that manifest themselves as different spatial distributions of 
the liquid in the connected pores. If the sets of macroscopic equilibrium states characterized by the 
distributions of filled and empty pores in the course of infiltration and defiltration were identical, 
the total thermal effect in the closed cycle would be equal to zero. In this case, the thermodynamic 
relationship (1T81) for the closed cycle is not satisfied. 

The closed cycle and, consequently, the complete transformation of the work into the heat was 
observed for a number of water- hydrophobized silica gel systems [7tfT7 t [T8 | l2T | l22] . In particular, the 
complete emptying of pores after the infiltration and the subsequent defiltration with a decrease in 
the excess pressure to zero was observed in [TTjlTH] for the KSK-G silica gel modified by n-alkylsilane 
molecules (n = 8, 16) grafted to the silica gel surface with a surface density higher than 2 nm -2 . 
Before the modification, the specific surface area, the pore volume, and the average pore radius 
for this silica gel were equal to 310 ± 20 m 2 /g, 0.95 cm 3 /g, and R = 5.2 nm, respectively. The 
values of these quantities after the modification are not presented in [TTPTB]. which complicates the 
analysis of the results obtained in these works. The closed cycle was also observed in [21] for the 
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water-C8W silica gel (Waters) system in which the silica gel was modified by n-alkylsilane with 
the chain length n — 8, the average pore radius R = 4.2, and the specific pore volume of 0.53 
m 2 /g in the temperature range 287 -j- 333 K. The authors established that the small nonoutflow 
(< 1%) takes place only at a temperature of 278 K. For the systems containing water and the 
Fluka 60 C8 silica gel, as well as the Zorbax Z4, Z8, Z18, and PEP10C18 silica gels, the closed 
cycle was observed in [7J|22]. In [7], the authors investigated the infiltration and defiltration in 
four porous media MCM41 with pores in the form of cylindrical channels. These media were also 
modified by n-alkylsilane with n = 8 and a surface density of 2.1 nm" 2 and had the average pore 
radii R = 1.3, 1.6, 2.0, 5.4 nm. The phenomenon of nonoutflow was observed only for the porous 
medium with R = 5.4 nm, whereas the other porous media with R = 1.3, 1.6, 2.0 nm underwent 
a closed infiltration-defiltration cycle. 

For the porous medium with specified parameters of the macroscopic structure of the pore 
space, the closed cycle with the complete defiltration, according to relationship (fT8|) . is possible 
only when the values of the quantities a, 5a, da/dT and d5a/dT fall in particular ranges. There- 
fore, relationship f fTB"]) requires a separate detailed quantitative analysis, which will be performed 
elsewhere. It should be noted that this analysis, in turn, necessitates the knowledge of the values of 
the quantities a, 5a, da/dT and d5a/dT; the porosity; the pore size distribution (which is changed 
after the modification |19|): and the quantities AQ and AA. Here, we restrict our consideration 
to the qualitative analysis of the available experimental data for which relationship (TLB"]) is not 
satisfied and the sum of the work and the heat in the cycle is equal to the change in the internal 
energy of the system. In this case, the mechanical work dissipated by the system is not equal to the 
total heat release, as it was described in |44| . The phenomenon of nonoutflow associated with the 
change in the internal energy is characteristic of the majority of the studied hydrophobic porous 
media and liquids, namely, water |7TO|[l^[2^[2^|2^[3^H35l[37] , aqueous solutions of salts [151122128], 
organic compounds, ethylene glycol [5J, alcohol [26J, and glycerol [45], as well as systems with the 
liquid metal [T|, Wood's alloy [I], and mercury [6Tjl6"2]. 

For water, the derivative involved in relationship (1181) has the value da/dT = 1.5 • 10~ 4 J/m 2 
K [53J. As the temperature changes from 293 K to 353 K, the surface energy a in accordance 
with this value of the derivative changes by 5%. Such a small change of the surface energy a 
in the system containing water and the modified Fluka 100 C8 silica gel appears to be sufficient 
that the closed cycle will transform into the cycle with a nonoutflow of more than 80% of water 
with a decrease in the temperature from 353 to 293 K [21]. This cycle is characterized by a change 
in the internal energy of the system (AE ^ 0). The quantity AE reflects both the reversible 
and irreversible changes in repeated infiltration-defiltration cycles. The reversible change can be 
associated with the adsorption of water on the modified surface |19[|20"]. This water evaporates 
already at room temperature (and more rapidly at an elevated temperature), and the first cycle 
with nonoutflow is reproduced 01211 EH]- The differential thermal analysis performed in [TJ] 
demonstrated that water evaporates at a temperature T = 373 K. The irreducible change of the 
internal energy AE is governed by the interaction of the liquid with the surface of the porous 
medium and depends on the maximum pressure during the infiltration [25J, the temperature and 
time of heating of the system [23], and the procedure used for preparing the surface after the 
modification [12] . The change of the internal energy AE also depends on the length (n) of the 
grafted n-alkylsilane molecule. It is worth noting that, for water and silica gels, an increase in n 
from 4 to 18 leads to a increase 5a, a decrease in the nonoutflow and a decrease in the change 
of the internal energy AE |7][8| [T&H20] . whereas for water and the modified medium, the closed 
cycle is observed for n = 1 and the complete nonoutflow takes place for n = 8 [32] . An increase in 
the concentration of the NaCl or CaCl2 salt in the aqueous solution in specific ranges leads to an 
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increase in the quantity 5a and a decrease in the change of the internal energy AE for the Fluka 
100 C8 silica gel |25lf2Tl 1281132] . A decrease in the concentration of ethylene glycol in the aqueous 
solution and the hydrophobic silica gel L23 results in the transition from the complete nonoutflow 
at a concentration c > 15% to the closed cycle (the nonoutflow is less than <5%) [5]. Unlike 
ethylene glycol, ethanol wets the modified surface of the Fluka 100 C8 silica gel. This brings about 
the adsorption of ethanol and an increase in the change of the internal energy |Ai?| |26| . 

In recent years, experimental data have been published on the infiltration of hydrophobic 
zeolites with water and aqueous solutions of salts [T0 | [TT | [3^14^0] . It has been revealed that, for 
systems containing water and silicalite 1 (OH), silicalite 1 (F~), and ZSM-5 zeolite, the hysteresis is 
not observed and the pressure dependences of the volume of the system measured for the infiltration 
and defiltration coincide with each other. This means that such systems exhibit properties of an 
elastic spring without dissipation at AE = in the cycle. With a cyclic change in the temperature 
from 358 to 318 K, these systems manifest a volume memory effect [TO] ITT]. For the water- 
MFI zeolite system, the dependences of the pressure on the volume and temperature during the 
infiltration remain unchanged when the rate of decrease in the volume changes by three orders 
of magnitude [38]. This implies that the Laplace- Washburn and Poiseuille equations are not 
applicable to the system under investigation. The infiltration-defiltration hysteresis was observed 
for a KC1 aqueous solution and the Zeolyst CBV-901 (HY) zeolite treated with SiCL [llj. This 
hysteresis depends on the nature of the anion (30] and, for the Y(ZY)-zeolite, on the nature of 
the cation (Li, Na, K, Cs); in this case, the infiltration pressure decreases with an increase in the 
cation radius [39] . 

In the analysis of experimental results obtained for zeolites, it is necessary to take into account 
that, for a channel (pore) diameter smaller than 1 nm, the liquid acquires properties of one- 
dimensional systems [38 | l64" j [65 ] . which differ qualitatively from the properties of the liquid in 
channels of larger sizes. This problem requires a separate analysis. 

In conclusion of the discussion of relationship (|18j) . we should note that it was derived under 
the assumption that 5R/R < 3 in the absence of correlations in the mutual arrangement of pores 
of different sizes and with the use of the model of pores as a system of randomly arranged spheres. 
Relationship (JT8]) for the closed cycle does not contain the average pore radius R; however, the 
value of R affects the porosity if and, hence, the number of neighbors z, the structure of the 
percolation cluster, and the quantities W{6) and W\{6). 

5 TEMPERATURE DEPENDENCES OF THE INFILTRA- 
TION AND DEFILTRATION PRESSURES 

Now, we analyze the signs of the derivatives da/dT and dSa/dT involved in relationship (TT81 
and the integrals of the quantities P(6) and + It follows from relation fllOp that the 

sign of the total thermal effect during infiltration is determined by signs of AQ p , AQ W values. The 
sign of the thermal effect due to the formation of the liquid-solid interface AQ p depends on the 
sign of d5cr/dT, which can be ascertained using known dependences of pressure at the beginning 
of infiltration p in and defiltration p out on temperature. The experiments carried out showed that 
for all investigated systems (modified silica gel - water, aqueous solutions of salts) the pressure at 
the beginning of infiltration changes not much as temperature increases dpi n /dT = 0, while the 
defiltration pressure increases as temperature rises [2~Tjl22 ,24j. Therefore, in order to determine the 
sign of the derivative d5a/dT, we calculate the pressures required for the infiltration of one pore 
on the shell of the infinite cluster and the defiltration of the liquid from an arbitrary pore under 
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the conditions SAi n = and 5A out = 0, which follow from relationship fl3j). By assuming that 
the quantities rj, R, and W do not depend on the temperature and taking into account that the 
infiltration and defiltration begin in the vicinity of the corresponding percolation threshold from 
relationships (fT3|) and (115]) . the derivatives dpi n /dT and dp out /dT can be written in the following 
form: 



dpin 


3 


dda 


dT 


~ R ' 


dT 


dpout 


3 


dSa 


dT 


" R ' 


dT 



Or 



[i-v)- 



It follows from eq. flE?} , that W{9 ~ 9 C ) > 0, Wi( 
it follows from eq. (1T9"]) . (ED]) that 



-1-6 
1-0 C ) >0, # c 



c/cr 



[i-v)~-^'vW(e~o c 



(19) 
(20) 

0.18. Since dp in /dT ^ 0, 
(21) 



Hence, the sign of dSa/dT is opposite to the sign of da/dT. The coefficient of the surface tension 
of the liquid-gas interface decreases as temperature increases, and vanishes at critical point so that 
da/dT < |66j. For water, value da/dT is —1.5 • 10 -4 J/m 2 K [66j. In this case, from eq. (12 Op find 
that for pressure at the beginning of infiltration dp out /dT > 0, which corresponds to experimental 
data [211123]. 

i 

Since the probability that a pore belongs to an infinite cluster is P(0) < 1 and J (6P(6) — 6)d6 < 

o 

0, it is necessary to analyze the dependences W(9) and W\{9) in order to determine the sign of 

the thermal effect in the course of the infiltration and defiltration. The behavior of the quantities 

l 

Wi(9) and W(9) and, hence, the sign of the integral J (Wi(9) + W(9))d9 depend substantially on 

o 

the porosity if. The figures shows the dependences W(9) and W\{9) calculated from relationships 
(112"]) and ( ]T6]) for different values of the porosity tp. 
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Figure 2: The dependences W(9) (solid line) and W\(l — 9) (dashed line) for porosity (p = 0.25 

It follows from the figures that the W{9) and W\{9) functions nontrivially depend on porosity. 
Since the filling of the porous medium begins only with the formation of the infinite cluster, 
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Figure 3: The dependences W(9) (solid line) and W±(l — 9) (dashed line) for porosity ip = 0.3 

































0.3- 
















0.2- 










/ 






0.1- 








/ 








0- 
-0.1- 
















\ 

s 


















0.2 


0.4 


0.6 

e 




0.8 





Figure 4: The dependences W(9) (solid line) and W\(l — 9) (dashed line) for porosity (p = 0.6 

W(9) = for 9 < 9 C = 0.18. At low porosity ip ~ 0.22 the structure of pores being filled is close 
to the fractal structure of the low-density infinite cluster near the percolation threshold, so that 
the growth of its surface is compensated by the decrease in the difference between the number of 
menisci in final and initial states of the filled pore at 9 < 0.3, since in concordance with §§§ there 
are few nearest neighbors of this pore z < 3 (Fig. [2]) . During further infiltration only the reduction 
in number of emerging menisci occurs (Fig. [2]), which reaches minimum at 9 w 0.68. If 9 = 1 the 
value W(l) is zero, since the filled porous medium meniscus is absent. 

For media with high porosity ip > 0.3 in which the number of nearest neighbours is z > 3, as 
the degree of filling 9 > 9 C increases, the infinite cluster of filled pores grows, accompanied by the 
growth of its surface and increase in the number of contacts of an empty pore with its neighbours 
on the shell of the infinite cluster filled with liquid. It leads to decrease in the difference between 
the number of menisci in final and initial states of the filled pore, so that as 9 increases the W(9) 
value reaches maximum and then decreases (Fig. OH]). In accordance with (fl2"|) it will continue 
until decrease in the difference between the number of menisci in final and initial states of the 
filled pore compensate the growth of area of the infinite cluster, which will lead to vanishing of 
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W(9) at 9 ~ 0.5 (Fig. [2J |3]). The further growth of the degree of filling leads to further decrease in 
the difference between the number of menisci in final and initial states of the filled pore while the 
growth of the surface of the infinite cluster of filled pores slows down. As a result W(9) reaches 
maximum at 9 ~ 0.7. The surface of the infinite cluster of filled pores will vanish at 9 = 1, which 
will lead to vanishing of W§B) (Fig. [HE]). The dependence of the W\(l — 6) value, which determines 
defiltration, is similar to the dependence of W(9), which determines infiltration. The difference in 
behavior of W(9) and W\(l — 9) is connected with the fact that liquid defiltration from the porous 
medium does not require the percolation cluster formation. Figures [2j|4] show that W(9 ~ 9 C ) > 0, 
Wx{9 ~ 1 - 9 C ) > 0. 

i 

The calculation of the integral J (Wi{9) + W(9))d9 shows that, with an increase in the porosity, 

o 

this integral increases from negative values at (p = 0.25 to positive values at (p = 0.6 and vanishes 
at if = 0.3. 

Expressions (TLB"]) . (fT5]) . (TL9]) . ([20]) and conditions 5A in = and 5A out = allow one to calcu- 
late the temperature dependences of the infiltration pressure Pi n (T) and the defiltration pressure 
Pout{T). Figure [5] shows the experimental data for silica libersorb 23 (silica gel KSK-G with the 
modification of 8-tier alkynesilan-C8), FlucalOO C8 [24], and S8W [21J filled by water and calcu- 
lated according to the considered correlation effects. Experimental data within the measurement 
error are described by linear dependences with different slopes. With increasing temperature, the 
pressure of infiltration decreases at about 10% and the defiltration pressure increases in the times 
for all the porous media. For these environments values of pressure p in and p out also differ at ini- 
tial temperature. In concordance with relationships 0131) and 0T5|) . dependences Pi n (T) and p ou t(T) 
are described by the quantities W(9 = 9$ > 9 C ), W\{9 = 1 — 9 ), and by dependences er(T) and 
5a (T). The dependence of <j(T) is known [66j. The value 5a in our experiments for libersorb 23 
was 5a = 22 ± 1 J/m 2 . The same value was taken for the other two porous media because of 
their similar modifications. The value d5a/dT was calculated from the condition ( l2~TT) . Quanti- 
ties W and W\ calculated according to expressions ([9]), ( fi~2l) . f lT6l) . The porosities estimated from 
the specific pore volumes according to the data taken from [2T1I24"] were as follows: tp = 0.33 for 
Libersorb 23, <p = 0.46 for Fluka 100, and <p = 0.53 for C8W. The quantities W{9 = 9 > 9 C ) and 
W\(9 = 1 — 9\, 9\ << 1) corresponded to the maximum values in the curves W(9) and W\(l — 9) 
for ip = 0.33, ip = 0.46, and <p = 0.53, respectively. 

As is also seen from Fig. El the calculated temperature dependences of the infiltration and 
defiltration pressures for the porous medium Libersorb 23 C8, Fluka 100 C8, and C8W infiltrated 
with water satisfactorily fit the experimental data. 

It follows from relationships ( JT3|) and ( !T5|) that, in the framework of the proposed model, the 
pressures of the beginning of the infiltration and defiltration p in and p out are proportional to 1/R if 
the quantities 77, W, and W\ are independent of R. This is consistent with the known experimental 
data |61I7]IT!? | |4"3" | I5T) [ I52"] . Such dependence is a consequence of the lack of correlations in the mutual 
arrangement of pores of different sizes in a model of randomly arranged spheres with a narrow pore 
size distribution 5R/R < 3, characteristic for silica gels. For the porous medium MSM 41 more 
than the strong dependence of the average radius [6l[7] due to the peculiarities of the fluctuation 
formation of the nucleus in a cylindrical channel (see (j3J)). 
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Figure 5: Dependences (a) pj n (T)and (b) p ou t{T) calculated for Libersorb 23 (curves 1), Fluka 
100 C8 (curves 2), and C8W (curves 3) infiltrated with water according to the calculations from 
relationships f|T3|) and (115)) . Points are the experimental data taken from [24"[|4"3] 

6 THERMAL EFFECT 

Expressions ([TO]) and ( TT2|) allow one to calculate the thermal effects observed during the infil- 
tration of a porous medium with a nonwetting liquid in different cases as functions of the porosity 
and surface energies of the liquid and the porous medium. 

It was found in [43j that when one-third volume of the porous medium from Sigma- Aldrich was 
filled by water, the temperature did not increase within the limits of error (< 0.1K). A maximum 
value of the elastic energy in performed tests was about 1 J. Authors |43] estimated the possible 
increase of temperature, in condition that the work (A) of filling is to raise the temperature. In 
experiments [43J at defined value A as 2.9 J, the temperature increase should be AT = 0.8K |43| . 
Figure [5] shows the calculated (from relationships ( 1T0|) and ( fT7|) ) dependences of the thermal effect 
on the degree of infiltration of the porous medium (in relative units) due to the formation of the 
liquid-porous medium surface AQ P (8)/A (Fig. [5b., lower curve) and menisci AQ W (8)/A (Fig. [5b., 
upper curve) and the corresponding dependence of the total effect AQ/A = (AQ P + AQ W )/A 
(Fig. [5b). In accordance with ( JTOj) . (fT2j) . the value of the thermal effect depends on the value of 
the integrals appearing in these relations, which, (see (191) . ([T2~]) and (fTB"]) ) depend on the porosity. 
Estimates show that in these experiments the porosity of silica gel after modification was ip 0.68. 
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The parameter T • d5a/dT was calculated from the condition (121j) . 
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Figure 6: Dependences of the thermal effect (in relative units) on the degree of infiltration 9 of 
the porous medium (</? f=s 0.68) with water: (a) the contributions due to the formation of the 
liquid-porous medium surface AQ P (9)/A (lower curve) and menisci AQ W (8) /A (upper curve) and 
(b) the total thermal effect AQ/A = (AQ P + AQ W )/A, A=2.9 J 

Note that W(9) is equal to zero at 9 < 9 C . Therefore, in Fig.[6]for comparison with experimental 
data, degree of infiltration is delayed on the horizontal axis, X shifted by an amount 9 C = 0.18 so 
that the value 9 = corresponds 9 = 9 C in the equation (TT2]) . 

It can be seen from Fig. E]that the thermal effect associated with the infiltration of the porous 
medium with the above parameters is small: at the maximum, it reaches ~ 0.15^4 for the degree of 
infiltration X ~ 0.17. For X « 0.27, the thermal effect vanishes. This is explained by the different 
origins of the contributions from the menisci and the pore surface to the total thermal effect, so 
that each contribution is one order of magnitude larger than the total thermal effect. 

Thus, in the performed experiments |UJ, upon the infiltration of the porous medium with 
water (the heat capacity is 4.2 J/gK [66]) for the heat release Q ~ 0.15v4 = 0.45 J, the maximum 
increase in the temperature is AT ~ 0.2K. In this case, it should be expected that, according to 
relationship (fTB"]) . the change of the internal energy of the system AE in the performed experiments 
upon the transition from the initial state to the final state during the infiltration and defiltration 
differs from zero and is comparable in the order of magnitude to the work expended for infiltrating 
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the porous medium: AE ~ A. This internal energy is equal to the difference between the energy 
of menisci in the pores on the shell of the infinite cluster of filled pores and the surface energy of 
the porous medium-liquid interface in the filled pores. The positive energy of system compression 
also give contribution to AE. 

Relationships ([9]) and (fl2|) make it possible to compare the values of the heat release observed in 
the experiments on the infiltration of the KSK-G (modification C16) porous medium with water [2j. 
According to the estimates, the medium porosity in the experiments performed in [2] is ip « 0.22. 
For this medium, we calculated the thermal effect associated with the infiltration. Figure [7J shows 
experimental data and the dependence of the thermal effect (in relative units) AQ{9)/ \Qq\ (where 
\Qq\ is the maximum heat release upon the infiltration of the porous medium with water; according 
to P|, Qo — 4 J/g) calculated from relationships fl9]) and (1120 . The elastic energy of compression 
of the porous medium and water is determined by water AQ U = —TAV ■ dK/dT, where AV is 
the specific decrease in the volume of the water-porous medium system and K = 0.166 mN/niK 
is the compressibility of water [2j. The maximum value of | AQ U \ is approximately equal to 6 J/g. 



Figure 7: Dependences of the heat release Q = AQ(8)/Qq on the degree of infiltration according 
to the measurements performed in [2] (points) and calculations from equations ([9]), 0120 with 
allowance for the compressibility of the porous medium with ip ~ 0.22 

Figure [7J also presents the components of the thermal effect measured in [2j : the dependence 
of the heat release associated with the compressibility of the system (curve 1) and the dependence 
of the heat release AQ p due to the change in the surface energy of pores (curve 2) and menisci 
AQ W (curve 3). It should be noted that the porous medium studied in [2] has a porosity <p ~ 0.22 
and, hence, its system of pores is located in the vicinity of the percolation threshold tp c = 0.18. 
Consequently, the infinite cluster of empty pores, which is required for the infiltration of the porous 
medium, is very sparse |10| and the probability (involved in relationship f fTUl) ) that a pore belongs 
to an infinite cluster is low P <C 1, therefore term contribution AQ P to the thermal effect is small. 
In this case, it follows from relationship f lTOj) that, in the infiltration of this medium, the decisive 
role will be played by the second term in expression (fTUj) . which containing integral from W(8) and 
corresponds to the contribution of menisci to the thermal effects (AQ W < 0). It is this circumstance 
that is responsible for the unusual thermodynamic properties of the system used in PJ. It can be 
seen from figure [7J that calculated dependences of total thermal effect is in good agreement with 
the experimental data. 
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Relationships (l9l)-ffT2l) make it possible to explain the thermal effect observed in the experiments 
[461155] on KSK-G (modification C16) with porosity (p ~ 0.4, on PEP100 (modification C18) and 
PEP300 (modification C18) with porosity (p ~0.65 and <p ~0.58 correspondingly. The average radii 
of the investigated porous medium are 6.5 nm, 5 nm and 15 nm for the KSK-G (C16), PEP100 and 
PEP300, respectively. The infiltration of the KSK-G silica gel was carried out at 308° K and the 
infiltration of the PEP100 and PEP300 were carried out at 298° K. Figures IHllTOl show experimental 
data [46,55j and dependences of the thermal effect (in relative units) calculated from relationships 
(ED-flHD for the KSK-G (C16), PEP100(C18) and PEP 300(C18). The filled volume was normalized 
and transformed to the degree of filling in calculations for dependences from relations (j9|)- (|T2|) . The 
elastic energy was estimated from the experimental data of dependence of pressure on the volume, 
taking into account the compressibility of the container and liquid. The parameter T • d5a/dT was 
calculated from the equation ([21]) as previously. 
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Figure 8: Dependence of the heat release on the degree of infiltration according to the measurement 
performed in [47] (points) for porous medium KSK-G and calculations from equations (p|)- (|T2|) with 
allowance for the compressibility of the porous medium ip ~ 0.4 

Thus, for the investigated porous media with tp = 0.4, (p ~ 0.58 and <p ~ 0.65 value of the 
thermal effect is negative and decreases with increasing the degree of filling, reaching a minimum 
at complete infiltration. It is indicates, that menisci give a significant contribution to thermal 
effects observed in [46j|47j. It can be seen from figures [7HTU1 that calculated dependences of total 
thermal effect is in good agreement with the experimental data. 

The measurements performed in [9UT9] during multi-cycle infiltration-defiltration process (num- 
ber of cycles ~ 1200) demonstrated that the increase in temperature in the systems containing 
hydrophobized silica gels and water per cycle is less than < 10 _2 K. The analysis of the experimen- 
tal data reported in [19] showed that, in the course of the cyclic infiltration-defiltration process, the 
hysteresis and, hence, the heat release in the first cycle are two times larger than those observed 
in the 1200th cycle. Thus, in the experiments conducted in [19], there occurs a partial nonoutflow 
of the liquid from the porous medium, which leads to changes in parameters and characteristics 
of the porous medium, such as the porosity, the average number of nearest neighbors, and the 
interfacial energy. In this case, it is necessary to perform a detailed analysis of the changes in the 
parameters of the porous medium in each cycle with inclusion of the thermal diffusivity of the 
porous medium and water in terms of the above-derived relationships, which is beyond the scope 
of our present work. 
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Figure 9: Dependence of the heat release on the degree of infiltration according to the measurement 
performed in |46j (points) for porous medium PEP300 C18 and calculations from equations (l9^- ffT2]) 
with allowance for the compressibility of the porous medium (p ~ 0.58 

7 CONCLUSIONS 

Thus, in this work, we have established the relation between the energy properties of the 
interface and the macroscopic properties of the pore space in a disordered porous medium and 
proposed the mechanism of energy absorption during the infiltration of the nanoporous medium 
with a nonwetting liquid. It turned out that thermal effects observed during the infiltration of the 
porous medium and the defiltration of the liquid from it can be positive and negative depending 
on the porosity and, hence, on the number of nearest neighbors of filled pores. 

The proposed mechanism is based on the inclusion of correlation effects during percolation 
infiltration of an infinite disordered porous medium. It is assumed that the infiltration of the 
porous medium is the result of growth of a percolation cluster consisting of filled pores through 
the attachment of empty pores (accessible to infiltration) to the shell of this cluster. Menisci 
appear and disappear in pores on the shell of the percolation cluster in the course of its growth, 
and these processes depend on the degree of infiltration. 

The above analysis is based on the representation of a system of pores in a porous medium in 
terms of the model of randomly arranged spheres. However, in this model, the correlations in the 
mutual arrangement of pores of different sizes are ignored and, hence, it is impossible to adequately 
describe the effects of blocking of the liquid in pores with large radii that are surrounded by pores 
with smaller radii. Therefore, in the framework of the model of randomly arranged spheres, we 
can describe the infiltration of disordered porous media only with a narrow distribution of pores 
over the sizes (5R) <C R. We assumed that pores in the porous medium have a regular (either 
spherical or cylindrical) shape, ignored the coordinate dependence of the surface energy of a pore, 
and operated actually with the average values of this energy on the surface of the corresponding 
pore. 

One of the consequences of the proposed mechanism is a condition that determines the class of 
systems for which a closed infiltration-defiltration cycle can exist. According to this condition, the 
initial and final states of the system coincide, the change of the internal energy is equal to zero, 
and the work expended for infiltrating the porous medium, which is determined by the area of the 
hysteresis loop, is equal to the thermal effect. Another consequence of the proposed approach is 
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Figure 10: Dependence of the heat release on the degree of infiltration according to the measure- 
ment performed in |46j (points) for porous medium PEP100 C18 and calculations from equations 
(|9])-( TT2]) with allowance for the compressibility of the porous medium ip ~ 0.65 

the dependence of the effective contact angle on the degree of infiltration of the porous medium 
with a liquid (see relationships (ITS]) and f)15p ). This dependence results from different paths on the 
way to the final state of the disordered porous medium during the infiltration (defiltration) with 
a liquid. 

Thus, the proposed approach makes it possible within a unified context to describe the tem- 
perature dependences of the infiltration and defiltration pressures for a porous medium with a 
disordered structure and the thermal effects associated with the absorption of the energy by "dis- 
ordered porous medium-nonwetting liquid" systems. 
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